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Abstract 

New data for the total cross section a{e~^e~ — > hadrons) in the charm and bottom 
threshold region are combined with an improved theoretical analysis, which includes 
recent four-loop calculations, to determine the short distance MS charm and bottom 
quark masses. A detailed discussion of the theoretical and experimental uncertain- 
ties is presented. The final result for the MS-masses, mc(3 GeV) = 0.986(13) GeV 
and mh{W GeV) = 3.609(25) GeV, can be translated into mdmc) = 1.286(13) GeV 
and ■mb{mh) = 4.164(25) GeV. This analysis is consistent with but significantly more 
precise than a similar previous study. 

PACS numbers: 12.38.-t 14.65.Dw 14.65.Fy 



1 Introduction 



The strong coupling constant and the quark masses are the fundamental input parameters 
of the theory of strong interaction. Quark masses are an essential input for the evaluation 
of weak decay rates of heavy mesons and for quarkonium spectroscopy. Decay rates and 
branching ratios of a light Higgs boson, suggested by electroweak precision measurements, 
depend critically on the masses of the charm and bottom quarks, rric and m?,. Last not 
least, confronting the predictions for these masses with experiment is an important task 
for all variants of Grand Unified Theories. To deduce the values in a consistent way from 
different experimental investigations and with utmost precision is thus a must for current 
phenomenology. 

Let us recapitulate the main ingredients of the sum rule approach which will be used in 
the following analysis. Originally the idea has been suggested for the analysis of the charm 
quark mass by the ITEP group [I] long ago. Subsequently the method has been developed 
further [2] and frequently apphed to the bottom quark. Most of these later analyses 
concentrated on using relatively high moments which are less sensitive to the continuum 
contribution and exhibit a very strong quark mass dependence [3l HI [5l [6l [7] . However, 
this approach requires the proper treatment of the threshold, in part the resummation 
of the higher order terms from the Coulombic binding and a definition of the quark 
mass adopted to this situation like the potential- or IS-mass [HI E], which subsequently 
has to be converted to the MS-mass. The low moments, in contrast, can be directly 
expressed by the short- distance mass at a scale around 3 GeV and 10 GeV for charm and 
bottom, respectively. The extrapolation to even higher scales, required for a number of 
applications is therefore insensitive to the larger corrections which would appear in the 
low energy region. 

A detailed analysis of rric and rrib, based on sum rules, has been performed several 
years ago [10] (see also Ref. [11]) and lead to mcirric) = 1.304(27) GeV and mhirri},) = 
4.191(51) GeV. This is still one of the most precise results presently available. During the 
past years new and more precise data for a{e~^e~ hadrons) have become available in the 
low energy region, in particular for the parameters of the charmonium and bottomonium 
resonances. Furthermore, the error in the strong coupling constant with its present valued 
as{Mz) = 0.1189 ± 0.0020 [12] (as compared to a,(Mz) = 0.118 ± 0.003 for the last 
analysis) has been reduced. Last not least, the vacuum polarization induced by massive 
quarks has recently been computed in four- loop approximation [T^ 115]; more precisely: 
its first derivative at = has been evaluated, which corresponds to the lowest moment 
of the familiar i?-ratio. A fresh look at the determination of the quark masses based on 
these new developments is thus appropriate and will be presented below. 

The outline of the paper is as follows: The basic assumptions of our approach are 
presented in Section [2l In Section [3] we discuss in detail the evaluation of the background 
and recall the theory input required for the quark mass determination. The measurements 
in the charm threshold region are discussed and applied to the determination of the charm 

^We adopt the central value from Ref. [H] but double the error such that the interval overlaps with 
the central value given by the PDG [13] . 
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quark mass in Section HI Similar considerations are used in Section [5] to obtain the bottom 
quark mass. Section [6] contains our conclusions. 



2 Generalities 

The extraction of mg from low moments of the cross section a{e^e^ — » QQ) exploits 
its sharp rise close to the threshold for open charm and bottom production and the 
important contributions from the narrow quarkonium resonances. At best the properties 
of the lowest bottomonium state, T(15'), can be evaluated in perturbative QCD. In general 
a differential description of the cross section close to threshold involves necessarily low 
scales, of O^agmq) or C(a^mg) and even non-perturbative contributions, arising, e.g., 
from condensates (G^), will enter. In contrast, by evaluating the moments 

- J^RQis), (1) 

with low values of n, the long distance contributions are averaged out and A4n involves 
short distance physics only, with a characteristic scale of order -^threshold = 2mg. Through 
dispersion relations the moments are directly related to derivatives of the vacuum polar- 
ization function at = 0, 



= ^ d? ""'^'^ 



(2) 



which can be evaluated in perturbative QCD. From dimensional considerations one obtains 
rriQ ~ {Ain)^ which implies 

rriQ 2n Mn 

for the relative error induced by the experimental uncertainties. (The theoretical and 
parametric errors will be described below.) Given 6R/R ~ 1.5 — 3%, a precision close 
to 10 MeV for rric and 20 MeV for rrify seems within reach. An analysis at this level 
obviously requires to control a variety of corrections, e.g., contributions from the non- 
perturbative condensate in the case of m^. Furthermore, a careful definition of the heavy 
quark production cross section per se is required. We will discuss this now in detail for 
charm production, with the generalization to bottom production being obvious. 

The narrow charmonium resonances J/"^, \&(2S') and the higher excitations will ob- 
viously contribute to the moments. Open charm production exhibits a sharp rise, nearly 
like a step function. Beyond the \E'(3770)-resonance a few oscillations are observed which 
quickly level out into a fairly flat energy dependence. Around and above approximately 
5 GeV the cross section is well approximated by perturbative QCD and mass terms can 
be considered as small corrections. The sensitivity to mg is, therefore, concentrated on 
the small region from J/\I' up to approximately 5 GeV. 
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Figure 1: Feynman diagrams contributing to R{s) at order a^. A secondary charm quark 
pair is produced through gluon sphtting. 

At first glance one might try to extract a{e~^e~ — >■ cc) by measuring the cross sec- 
tion with tagged charmed hadrons, thus ehminating the contribution from hght quarks. 
However, such data are presently only available for selected narrow energy regions. Even 
more important, this approach is also not tenable from the conceptual side. A sizeable 
fraction of ^(3770) decays (not to speak of J/ip or ip') proceeds to non-charmed final 
states, nevertheless its total contribution should be attributed to the moments. On the 
other hand, processes with "secondary" charm production arising from gluon splitting (cf. 
Fig. [11(a)) should rather be assigned to the light quark cross section. Also cuts through 
singlet diagrams (cf. Fig. [2]) which first appear in order and where the fermion lines 
represent light or heavy quarks, exhibit thresholds at = 0, 4m^ and 16m^, and these 
contributions cannot be attributed to a specific quark flavour. For the determination of 
mg we therefore rely on the measurement of the total cross section, which is calculated 
through the absorptive parts of vacuum polarization diagrams. Various contributions to 
the cross section, arising from different cuts through the same diagrams, will always be 
treated together. For the moments Ain we only consider those contributions to n(g^) 
which are analytic for < 4mQ and exhibit a cut starting at = 4mQ. They result 
from diagrams with one massive quark loop connecting both photon vertices and will 
be denoted by IIq and Rq, respectively. The remaining diagrams with cuts starting at 

= will be treated as "background" . This includes all diagrams with one light quark 
loop connecting both photon vertices (including those with internal massive quark loops) 
and all singlet diagrams. It will be shown below that the background defined in this way 
exhibits a fairly flat behaviour across the threshold, despite the fact that individual cuts 
with branching points at = 4mQ exist. 

In the remaining part of this Section we fix the notation and give precise definitions of 
the quantities used in the paper. It is convenient to normalize the radiatively corrected 
hadronic cross section relative to the point cross section and to define the ratio 

^, , cr(e+e^ —i> hadrons) 

Ris) = ^ (4) 

where apt = 47ra^/(3s). As an inclusive quantity R{s) is conveniently obtained via the 
optical theorem from the imaginary part of the polarization function of two vector currents 
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Figure 2: Sample singlet diagram contributing to n(g^) which arises at order for the 
first time. The cuts (indicated by the dashed lines; not all possible cuts are shown) 
represent contributions to R{s). 

via 

R{s) = 127r Im [n(g2 = s + ie)] , (5) 

where Il{q^) is defined through 

{-q'9,. + q,q.)W) = ^ / dx e^«^(0|rj^(x)jt(0)|0) , (6) 

with being the electromagnetic current. 

The perturbative expansion of R{s) can be written as 



Ri^) = EE(^)"^f^ + ^sing, (7) 



where the summation is performed over all active quark flavours i. (The small singlet 
contribution starts to contribute at order a^.) For a comprehensive compilation of the 
individual pieces we refer to [HI [T71 118] where also explicit results are given. The full 
quark mass dependence is available up to order [20] . For the non-singlet term R^^^ 

the first three terms of the high-energy expansion are known [21] . 

The predictions for R{s) which are used to calculate the background and the continuum 
term far above charm or bottom threshold, respectively, are based on Eq. ([7]) where the 
up, down and strange quark masses are taken as massless. For the charm and bottom 
quark the respective pole masses are chosen as input. If not stated otherwise we will use 
the following parameters for the evaluation of the background 

af\Mz) = 0.1189 ±0.002, 

Mc = (1.65 ±0.15) GeV, 

Mb = (4.75 ±0.20) GeV, (8) 
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which cover the full range of all currently accepted results. In Eq. ([8]) and refer 
to the pole masses. 

At several places of our analysis the renormalization group functions and the matching 
conditions for the masses and the strong coupling are needed in order to get relations 
between different energy scales. The corresponding calculations are performed using the 
package RunDec |22j . 

3 The background 

As stated above, we distinguish three energy regions: First, the region of the narrow 
resonances J /if) and ip' , second, the "charm threshold region" starting from the D-meson 
threshold at 3.73 GeV up to approximately 5 GeV, where the cross section exhibits rapid 
variations and, third, the continuum region where pQCD and local duality are expected 
to give reliable predictions. In this last region the cross section is mainly sensitive to 
and fairly insensitive to rric- 

For the present analysi^ we use the data from the BES collaboration published in 
2001 [23] and 2006 [21]. Whereas the older data cover the whole range from 2.0 GeV 
to 4.8 GeV the newer ones provide a very precise measurement in the region around 
^/s 3770 GeV. In Fig.[3]the BES-results are shown together with the measurements from 
the MD-1 [25] and CLEO [26] experiments between 4.8 GeV and 10.52 GeV. As is evident 
from Fig. [3l pQCD with af\Mz) = 0.1189 provides an excellent description of all recent 
results. For example, the recent BES-value of i^'^^P = 2.141±0.025(stat.)±0.085(syst.), de- 
termined in the interval [3.650GeV, 3.732GeV] is in excellent agreement with the theoreti- 
cal prediction of i?(3.700 GeV) = 2.161±0.017 and the same is true for i?^^P(10.52 GeV) = 
3.56 ± 0.01 ± 0.07 from CLEO [25] to be compared with i?(10.52 GeV) = 3.548 ± 0.012. 

Up to C(a^), the prediction for charm quark production incorporates the full quark 
mass dependence. Starting from order also the charm quark mass dependence of 
"secondary" charm production has to be taken into account. This includes diagrams 
of the type in Fig. [T](a) as well as those from Fig. [D^b). In addition we include 0{a^) 
terms from the expansion in (M^/s)"" with n = 0, 1 and 2 [TOl I2T] . Last but not least 
contributions from virtual c quarks {^/s < 3.73 GeV) and b quarks (-^/i < 10.52 GeV) 
are included. Since our formulae are expressed in terms of a[^\ the latter are suppressed 
~ (as/7r)^s/(4M^) and decouple for s -C 4M^. These are included in the forth column of 
Tab. [H Pure QED final state radiation is tiny and taken into account for completeness. 
(For a related analysis at = 10.5 GeV see [27j.) 

Let us now discuss in detail the evaluation of the "background" . The results will be 
used to extrapolate Ruds from the region below charm threshold up to 4.8 GeV. Here 
we provide all formulae relevant for the charm threshold region; the modifications to the 

2 We limit this analysis to the results from BES [SSlEl], those from MD-1 [25] and from CLEO [2^ with 
systematic errors of typically 4.3% (BES [23]), 4% (MD-1 [25|) and 2% (CLEO [26]). Older measurements, 
in particular those from SPEAR and DORIS, are consistent with the new results. However, with their 
significantly larger errors they do not provide additional information. 
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Figure 3: -R(s) for different energy intervals around the charm threshold region. The 
solid line corresponds to the theoretical prediction, the uncertainties obtained from the 
variation of the input parameters and of /i are indicated by the dashed curves. The inner 
and outer error bars give the statistical and systematical uncertainty, respectively. 



bottom case are obvious. 

Below 3.73 GeV only u, d and s quarks are produced. To allow for a smooth transition 
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through the threshold the effective number of flavours will be chosen to be n/ = 4 and 
virtual charm quark effects are taken into account (for a compilation of the relevant 
formulae see Refs. O [El [H] ) . 
Specifically we define 



-^background — Ruds ~l~ -Ruds(cb) ~^ -Rsing ~l~ RqED 

where we take for the purely light degrees of freedom 



(9) 



uds 



3 E Ql 

i=u,d,s 



1 



71 



1.640 - 2.250L 



(-10.28 - 11.38L,^ + 5.063L 



SfJ, J 



a. 



n 



(10) 



with = ln(s//i^) and = cti^''(A^)- Contributions with virtual or secondary charmed 
quarks appear first in order and in this order they are known in closed analytical 
form [SHlinilSO]: 



R 



uds(cb) 



'^-"'uds(c) 



i=u,d,s \ ^^ J 



6R 



(2) 

uds(b) 



3 E Qi y{M^,s) 

i=u,d,s 



where and can be found in Eqs. (E.15) and (E.16) of Ref. [18] o The first three 
expansion terms for small and large center-of-mass energy are given by 



3<4M| 



3 E 

i=u,d,s 



+ 



[mi) 



SR. 



(2) 

uds(c) 



3 E Q' 

i=u,d,s 



l (L,^ - L,mJ + (0.06519 - 0.01481L,M. 
6 M^ 



-0.001231 + 0.0003968L,mJ + • • • 

1 /M'^V 
0.1153 + -L,^ + I V ) ("^-^"^^^ + + 



(12) 



(13) 



5i?uds(c) exhibits a logarithmic divergence for small s, which could be removed by matching 
to a theory with three effective flavours. However, restricting ourselves to y/s > 2 GeV 
this logarithm is numerically not large and we can restrict ourselves to Uf = 4. The 
corresponding a^-terms with one or two internal heavy quark loops, which in Eq. flTTl) 



■^Note that in Ref. [18] and are combined such that n/ = 3 for s < 4Af^ and Uf = A for s > 4M^ 
whereas we always have n/ = 4. 
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(3) 

is parameterized via ^-R^dsCc) only known in approximative form in the regions s <^ 



and s > [21jfl 



''-"'uds{c) 



s<4M2 



6K 



(3) 

uds{c) 



s>4M2 



3 E 



0.1528 + 1.005L,„ - 0.7222L:,, + 0.6944L,„L, 



- 1.005L,m. + 0.02778LL/, + ^ 



-0.1568 -0.1709L 



+0.03210L,^L,M. + 0.1996L,M. - 0.03724L,V, 
+ [^"j (-0.007178 + 0.0009988L,^ - 0.00006614L,^L,m, 



-0.0002714L,M, + 0.0003830LtM + • • • 



3 E 

i=u,d,s 



(14) 



+ 



M2 



- 1.236 + 1.819L,^ - 0.7222L^^ + (-6.476) 
10.58 + 1.979L,^ + 17.91L,M^ - 2L2^^ 



4.167L,^L,A/J + ... 



(15) 



with LsM^ = ln(s/M2). 

In Fig. m the exact result for 5-Ruds(c) is compared with the approximations of Eqs. f|T2l) 
and f|T3l) . As can be seen, the approximations work excellent even fairly close to threshold, 

/ON 

above and below. The corresponding plot for 5-Ruds(c) shown in Fig. [51 The dash-dotted 
curve corresponds to the approximation in Eq. (fl5l) and the dashed one to Eq. (fT^ . 
One can see a good agreement of the two approximations for y/s^3.5 GeV. Since for 
y/s^3.5 GeV the lower dashed curve should be very close to the exact result we adopt 
the result in Eq. f|T^ for the whole interval 2 GeV<y/s<5 GeV. It is interesting to note 



that the numerical effect of 6R 



(3) 



?(2) 



uds(c) comparable to the one of 5-Ruds(c)- 



,. However, their 



overall size is small, as can be seen in Tabs. [T] and [2] (see below) . 

The singlet contributions start at order and are fairly small. In the limit s 3> 4M^ 
-Rsing can be cast into the form [18] 



R, 



sing 



5>4M2 



3(^ 

TT 



( E Qi] (-0.4132) +Q, E 

\i=u,d,s,c j i=u,d,s,c \ / 



(16) 



There are no quadratic mass terms in the singlet contribution. Both terms in Eq. (|T6 
are proportional to Ql, since the sum over the light quark charges happens to vanish. 



*The corresponding result can easily be extracted from Eqs. (E.22), (E.25) and (E.26) of Ref. [T8] . 
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Figure 4: SR^^J^^^^-^ for 2 GeV < ^/s < 5 GeV. Both the exact result (sohd hne) and the 
approximations including terms of order [s/MlY (dashed) and {M^/sY (dash-dotted), 
respectively, are shown and n"^ = s has been chosen. 

The expansion for small ^/s has been considered in Refs. [311 E2] in the context of 
the Z boson decay. It is straightforward to adopt Eq. (20) of [31] to our situation for 
i?sing. It turns out that the lowest four expansion terms for s -C 4M^ are proportional to 
Qu + Qd + Qs = 0- Thus the first non-zero contribution is proportional to Ql and arises 
at order {s/M^y. This term can be found in Refs. 021 [32] • However, its contribution is 
numerically very small and can safely be neglected in our context. 

Combining the information about the small and large energy region for -Rging one 
obtains that its numerical contribution is between — 0.55(as/7r)^ and zero and can thus 
be safely neglected. Similar arguments also hold for the bottom case. 

The magnitude of the most important terms contributing to -Rbackground is shown for 
several characteristic energies in Tabs. [Hand [2] for both the charm and the bottom region. 
In Tab. [3] they are compared to a few selected experimental results. Note that in contrast 
to -Rbackground & three-(four-)flavour theory has been used for the y/s values below (above) 
the charm threshold region. All theory predictions are given for the reference values 
specified in Eq. ([S]). 
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Figure 5: SR^^J^^^^-^ for 2 GeV < y/s < 5 GeV. The approximations including terms of 
order {s/M^Y (dashed) and {Ml/ s)"^ (dash-dotted), respectively, are shown and 1x^ = 3 
has been chosen. 

The band in Fig. [3] and the uncertainty in Tab. [3] is obtained from the parametric 
errors and the theoretical uncertainty, the latter obtained by varying the renormalization 
scale between /x^ = s/4 and /i^ = 4s. These were added quadratically. The excellent 
agreement between prediction and measurement in the continuum justifies to use the 
background formula, normalized to data below threshold also in the region above to 
obtain as remainder R^. which will be used for the quark mass determination. 

In passing let us mention that the formulae in combination with the R measurement 
between 3 GeV and 10.5 GeV can be used to determine a fairly precise value for a^, as 
demonstrated in Ref. p^ . 
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Table 1: Predictions for -Rbackground in the charm region for selected energies, separated 
into various pieces as defined in the text. For completeness we also give the values of 
aW(v/i). 



(GeV) 




-Rudsc 


-Rudsc(b) 


-^background 


4.80 


0.2186 


3.7617 


0.0010 


3.7649 


5.50 


0.2102 


3.6828 


0.0003 


3.6852 


6.20 


0.2034 


3.6378 


-0.0003 


3.6396 


6.90 


0.1977 


3.6092 


-0.0006 


3.6107 


7.60 


0.1928 


3.5896 


-0.0009 


3.5908 


8.30 


0.1886 


3.5753 


-0.0011 


3.5764 


9.00 


0.1849 


3.5644 


-0.0012 


3.5654 


9.70 


0.1816 


3.5559 


-0.0013 


3.5567 


10.40 
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3.5489 


-0.0014 


3.5497 


11.10 


0.1759 


3.5431 


-0.0014 


3.5438 


11.80 


0.1735 


3.5382 


-0.0014 


3.5389 



Table 2: Predictions for -Rbackground in the bottom region for selected energies, separated 
into various pieces as defined in the text. For completeness we also give the values of 

«f(v^)- 



4 Charm quark mass determination from the thresh- 
old region 



Our determination of the charmed quark mass follows closely Ref. [10]. It is based on 
the direct comparison of the theoretical and experimental evaluations of the contributions 
to the derivatives of the polarization function Ilc{q^) (as defined in Eq. ([6])), the former 
evaluated in perturbative QCD, the latter through moments of the measured cross section 
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(GeV) 


2.00 


3.65 


3.732 


4.80 


9.00 


10.52 




2.209(91) 


2.161(18) 


2.160(17) 


3.764(64) 


3.564(17) 


3.548(12) 




2.18(7)(18) 


2.157(35)(86) 


2.156(86)(86) 


3.66(14)(19) 


3.62(7)(14) 


3.56(1)(7) 


Experiment 


BES [23] 


BES [Ml 


BES [Ml 


BES 23J 


MD-1 


CLEO [26 



Table 3: Comparison of the theory predictions for R{s) with the experimental results at 
a few selected values for i/i. 



for charm production in electron-positron annihilation. In its domain of analyticity Ilc{q^) 
can be cast into the form 

nc(g') = QlT^,T.Cnz\ (17) 

16vr^ „>o 

with Qc = 2/3 and z = /{Ami) where rric = nicdj?) is the MS charm quark mass at the 
scale /i. The perturbative series for the coefficients Cn in order was evaluated up to 
n = 8 in Ref. [HI [20] (recently [31] even up to n = 30). The four- loop contributions to Cq 
and Ci were calculated in Ref. [HJ [15] . The coefficients depend on and on the charm 
quark mass through logarithms of the form /^^ = ln(m^(/i)//i^). They can be written as 




In Tab. 1 the coefficients C^) are given in numerical form up to n = 4. The terms C^^^ 
with n > 2 are still unknown. To estimate the theoretical uncertainly from this source 
we will use the following bounds on C^^'^^: 

- 6.0 < °^ < 7.0 , -6.0 < C-f °^ < 5.2 , -6.0 < °^ < 3.1 . (19) 

The lower limit C^^^ > —6.0 is based on the observation that for fixed order in the 
absolute values of C^^\ C^^^ and C^^ in only one case exhibit a slight increase, in general, 
however, they tend to decrease with increasing n > 1 quite markedly. The upper limit 
is estimated from the dependence on the order in with n fixed, where we assume at 
maximum a geometric increase and thus C^^°^ < {C^'^^)'^ / C^^\ In our numerical analysis 
we also include the two-loop QED corrections to Cn which can easily be obtained from 
C^^\ It leads to a negative shift in rric by less than 2 MeV. 

Tab. m and Eq. ([T9|) essentially constitutes our theoretical input. We define the mo- 
ments 

127r2/d\" ,2, 



(20) 

r2=n 



12 



n 


MO) 








M-21) 


^.(22) 


Mm 


/^(31) 


Mm 


M33) 


1 


1.0667 


2.5547 


2.1333 


2.4967 


3.3130 


-0.0889 


-5.6404 


4.0669 


0.9590 


0.0642 


2 


0.4571 


1.1096 


1.8286 


2.7770 


5.1489 


1.7524 




6.7216 


6.4916 


-0.0974 


3 


0.2709 


0.5194 


1.6254 


1.6388 


4.7207 


3.1831 




7.5736 


13.1654 


1.9452 


4 


0.1847 


0.2031 


1.4776 


0.7956 


3.6440 


4.3713 




4.9487 


17.4612 


5.5856 



Table 4: Coefficients of the photon polarization function in the MS scheme as defined in 
Eqs. ( fT7|l and (fTSjl . = 4 has been adopted which is appropriate for the charm threshold. 



which leads to 

As demonstrated in Ref. [10], high moments (n > 4) are less suited for the mass extraction. 
The analysis will therefore be limited to n < 4 and the results for n = 4 will only be 
presented for illustration. 

With the help of the dispersion relation we establish the connection between the 
polarization function and the experimentally accessible cross section Rc{s). In the MS 
scheme 

^ciq') = S-J^^^^.+Ql-^A^ (22) 
ll'K'^ J s[s — q^) lovr^ 

which allows to determine the experimental moments in Eq. ([1]). Note, that the last term 
in Eq. fl22l) which defines the renormalization scheme disappears after taking derivatives 
with respect to q^. Equating Eqs. fl?Il) and ([T]) leads to an expression from which the 
charm quark mass can be obtained: 

l/(2n) 



"^c(/i) = ^ ( ) • (23) 



Non-perturbative contributions to the moments have been evaluated in [U |35] and, 
in view of their smallnes^, were neglected in Ref. [10]. With the further reduction of 
theoretical and experimental uncertainties this approximation has to be reconsidered. 

The non-perturbative contribution from gluon condensate is given by [H [35] 

^ (^(^G'^) .„(!-. ^S..). (24) 

with 

2n + 2r(4 + n)r(|) 

15 r(4)r(|+n 



(2n + 4)(^-/^J , (25) 



''This observation was confirmed in Ref. 1361 
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^(2^) 


Mvt(GeV) 
ree(keV) 
(a/a(M^))2 


3.096916(11) 
5.55(14) 
0.957785 


3.686093(34) 
2.48(6) 
0.95554 



Table 5: Masses and electronic widths [13] of the narrow charmonium resonances and 
effective electromagnetic coupling [38] at the appropriate scales. 



and bn taken from Tab. 2 of Ref. [35] : 

135779 , 1969 , 546421 , 661687433 , , 

bi = , &2 = , &3 = , &4 = • (26) 

12960 168 42525 47628000 ^ ^ 

For the numerical value of the condensate we take [37J 

— = 0.006 ± 0.012 (GeV)^ (27) 

TT / 

For the evaluation of 6J^^ we use as{3 GeV) = 0.254, fi = 3 GeV and nic = 1.3 GeV. 
The result for 6Ai^ are listed in Tab. [6] and will be subtracted from 6Ai^^ which enters 
Eq. (1231) . Their impact on rric is still negligible, although the contribution to the error 
starts to become relevant for n > 3. 

Let us now turn to the three different contributions which enter the right-hand side of 
Eq. ([1]): the J/\E' and \Ef' resonances (A^J^f''), the charm threshold region between 2M£)q ^ 
3.73 GeV and ^il = 4.8 GeV as measured by the BES experiment [23l[2l] {M^^'^''^), and 
the continuum contribution above si (A^™°*). 

The two narrow resonances are treated in the narrow width approximation 

R'^'^is) = Q^^^^e Sis -Ml). (28) 




where we use the parameters listed in Tab. [51 

In the charm threshold region (which includes \Ef(3770)) we have to identify the con- 
tribution from the charm quark, i.e. we have to subtract the parts arising from the light 
u, d and s quark. (As discussed in Section [31 charm production through secondary gluon 
splitting is treated as part of the background, the same applies to the singlet contribu- 
tion.) Technically this is done by determining a mean value for -Rbackground (cf. Eq. ([9|)) 
from the comparison of theoretical predictions and the BES data between 2 GeV and 
3.73 GeV and using the theoretically predicted energy dependence to extrapolate into the 
region between 3.73 GeV and 4.8 GeV [39]. The mean value at the energy just below the 
i5Z)-threshold is given by 

-^background = -R*^(s_) ri_ , (29) 
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with 

n. 



1 



where the sums run over the s-values where experimental data, R^^^^ are available. In 
our case the upper bound is given by ^/sZ = 3.73 GeV and R^^ corresponds to the the- 
oretical predictions from -Rbackground(s). In a next step we multiply -Rbackground(s) with 
the normalization factor n_ and subtract the result for each s-value in the energy region 
[3.73 GeV, 4.8 GeV] from the data before the integration is performed. The final result for 
the charm quark mass mc(3 GeV) changes by about —5 MeV if we do not take this renor- 
malization into account. The effect of the energy dependence of -Rbackground has practically 
no effect on the final result (less than 1 MeV). The normalization factors n_ = 1.038 for 
the 2001 data [23j and n_ = 0.991 for the 2006 data [21] can be considered as systematic 
offset of the experimental R measurement, as determined from the measurement below 
threshold and is determined under the assumption that the energy dependence of the true 
i?-ratio is indeed given by R^^. The magnitude of this factor is consistent with the quoted 
systematic error of 4.3% [23] and of 4.0% (4.9%) outside (within) the ^(3770)-region I24j 
which is kept in the subsequent analysis. The statistical error of n_ is negligible. 

In the continuum region above y/s = 4.8 GeV data are sparse and imprecise. On the 
other hand, pQCD provides reliable predictions for R{s), which is essentially due to the 
knowledge of the complete mass dependence up to order [19] and the dominant terms 
of order [21] . Thus in this region we will replace data by the theoretical prediction for 
R{s) as discussed in the beginning of Section [31 As emphasized above, the cross section 
is only weakly mc-dependent in this region. 

In Tab. [HI we present the results for the moments separated according to the three 
different contributions discussed above. The error of the resonance contribution is due to 
the uncertainties of the input parameters. The error of the charm threshold contribution 
is dominated by the correlated normalization error of approximately 4.0% to 4.9% of the 
BES data. 

To estimate the error on A^™"^* we varied the input parameters as stated in Eq. ([8]) and 
the renormalization scale between ■\/i/2 and 2^/s. The errors of the three contributions 
are added quadratically. It is instructive to compare the composition of the experimental 
error for the different moments. Generally speaking, it is dominated by the resonance 
contribution, specifically by the 2.5% uncertainty in the leptonic widths of the J/'^ and 
\E''. Compared to the analysis of [TU] the results are consistent. However, the present error 
amounts to less than half of the previous one. Already for the previous analysis the im- 
provement in the cross section measurement due to BES (from about 10-20 % systematic 
error down to 4.3%) was important. This is particularly true for the moment n = 1. The 
recent data for the \E'(3770)-region have lead to further improvement. The parametric un- 
certainties (from as and Mc) and the residual /^-dependence which affect Al™'^* are small. 



E 
E 



a 



1 



2 ' 



(30) 



15 



n 


^res 
'■n 


thresh 


A^cont 


A^cxp 






xlo("-i) 


xl0("-i) 


xlo("-i) 


xlo("-i) 


xlo("-i) 


1 


0.1201(25) 


0.0318(15) 


0.0646(11) 


0.2166(31) 


-0.0001(2) 


2 


0.1176(25) 


0.0178(8) 


0.0144(3) 


0.1497(27) 


0.0000(0) 


3 


0.1169(26) 


0.0101(5) 


0.0042(1) 


0.1312(27) 


0.0007(14) 


4 


0.1177(27) 


0.0058(3) 


0.0014(0) 


0.1249(27) 


0.0027(54) 



Table 6: Experimental moments in (GeV)~^" as defined in Eq. ([T]) separated according 
to the contributions from the narrow resonances, the charm threshold region and the 
continuum region above a/s = 4.8 GeV. In the last column the contribution from the 
gluon condensate is shown. 



n 


me (3 GeV) 


exp 


as 


/i 


np 


total 




mc(mc) 


1 


0.986 


0.009 


0.009 


0.002 


0.001 


0.013 




1.286 


2 


0.979 


0.006 


0.014 


0.005 


0.000 


0.016 


0.006 


1.280 


3 


0.982 


0.005 


0.014 


0.007 


0.002 


0.016 


0.010 


1.282 


4 


1.012 


0.003 


0.008 


0.030 


0.007 


0.032 


0.016 


1.309 



Table 7: Results for mc(3 GeV) in GeV obtained from Eq. fl23l) . The errors are from 
experiment, a^, variation of /x and the condensate. The error from the yet unknown 
four-loop term is kept separate. Last column: central values for mdmc)- 



The higher moments (in fact already for n above two) are increasingly dominated by the 
resonance contributions with their 2.5% uncertainty. The large uncertainty in {ag/ixG'^) 
starts to matter for n > 3. 

We use the results of Tab. Altogether with Eqs. fl2^ and ^IM in order to obtain in a first 
step mc(3 GeV). The results are listed in Tab. [7] together with the total uncertainty and 
its decomposition into the individual contributions from the experimental moments, the 
variation of as (cf. Eq. ([8])), the renormalization scale /i = (3±1) GeV and the condensate 
which are all added quadratically. To determine the /i- dependence, we evaluate Eq. (1231) 
at /i = 2 GeV and 4 GeV. Subsequently, using the renormalization group equation to 
four-loop accuracy, we evolve the result for mc(/i) to the reference point = 3 GeV and 
compare with the original value. 

The first moment is available in four-loop approximation, and the corresponding /i- 
dependence of the result is completely negligible. It is interesting to anticipate the corre- 
sponding /x-dependence for the four-loop result also for the higher moments. The corre- 
sponding results are displayed in Tab. [71 In the eighth column we show the effect of the 
unknown coefficients C"^^"-* {n > 2) where assumptions were made on their respective nu- 
merical growth as discussed above. For the extraction of mc(3 GeV) we assume C^^^ = 0. 
For SC^^^ we take the larger of the deviations estimated in Eq. (fT9|) . 

Let us compare the results displayed in Tab. [71 with those from [10]. For the moment 
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with n = 1 the main improvement originates from the new experimental results for the 
electronic widths of the J/\Ef and \Ef'. Some additional improvement is due to the new BES- 
data [24]. The four-loop term leads to a further reduction of the theoretical uncertainty 
— in view of the dominant experimental uncertainties it only leads to consolidation of the 
theoretical prediction. The error composition in the analysis based on = 2 is similar, 
as far as the experimental input is concerned. However, in this case the error from the 
i?-measurements is somewhat smaller, the error induced by the uncertainty in is larger. 
Furthermore, the error from 6C2^^^ has to be taken into account. Although at first glance 
this error looks small, a check of the assumptions on C2^^^ would be extremely useful. 

The analysis based on n = 3 leads to a result consistent with the one for n = 1 
and n = 2. In this case, however, the theoretical uncertainty estimated through the n 
dependence starts to become important and the total error is therefore larger. 

For the subtraction of the light-quark continuum we have assumed pQCD to be exact. 
Deviations from this idealized situation have been estimated in Ref. [40j. Based on the- 
oretical and phenomenological considerations it is assumed, that the oscillations of R{s) 
around the perturbative result, which are observed in the low-energy region and reflect the 
low-lying resonances, are damped by a power-law or an exponential factor. The following 
two correction factors were suggested in [ID] 

/i = 1 + 1.22 sin (2pyi - 6) , (31) 
with p = 3 GeV~^, 5 = 1.32 and, alternatively, 

f, = 1 -1.24exp(-^) sin 3.08) , (32) 

with 0"^ = 2 GeV^, B = 0.5 and Nc = 3. Potential effects from this modification are 
estimated by including the factors /i or /2 into the analysis]^ In the case of /i we also 
vary the phase 6 between and vr. 

The change in the contribution from the region 3.73 GeV < ^/s < 4.8 GeV is about 
1.2% for fi and n = 1 and significantly smaller in the case of f2- The 1.2% translate to a 
shift of —1 MeV for mc(3 GeV) which is completely negligible in the present context. 

The results of Tab. [7] are also shown in Fig. [6] where the central value and the uncer- 
tainties of mc(3 GeV) are plotted for = 1, 2, 3 and 4. For each n the results are shown 
for one- to four-loop theory input. For the 0{al) term the exact result is used for n = 1 
whereas for n = 2, 3 and 4 we use C^^^ = and the error estimates discussed above. It 
is nicely seen that the results for mc(3 GeV) further stabilize when going from three to 
four loops. At the same time the uncertainty is considerably reduced. Furthermore, the 
preference for the first three moments is clearly visible. Also the analysis for n = 2 and 
n = 3 leads to small errors, even if we include to contribution from the yet uncalculated 
dC^'^\ We emphasize the remarkable consistency between the three results which we 
consider as additional confirmation of our approach. 

^We thank A. Vainshtein for raising this question. 
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Figure 6: mc(3 GeV) for n = 1,2,3 and 4. For each value of n the results from left 
to right correspond the inclusion of terms of order a^, a], and in the coefficients 



Cn (cf. Eq. ( ITSl) ). Note, that for n = 3 and n = A the central values and uncertainties 
can not be determined with the help of Eq. ( |23l) in those cases where only the two-loop 
corrections of order are included into the coefficients Cn as the equation cannot be 
solved for mc(3 GeV). 



The result based on the moment with n = 1 is evidently least sensitive to non- 
perturbative contributions from condensates, to the Coulombic higher order effects, the 
variation of fi and the parametric dependence. The results for n = 2 and n = 3 are 
practically identical. Since the error in this case is somewhat larger, and furthermore 
affected by the estimate for SC^^\ we adopt 

mc(3 GeV) = 0.986(13) GeV, (33) 

as our final result. Transforming this to the scale- invariant mass mc(mc) [22] including 
the four-loop coefficients of the renormalization group functions one finds 

mc(mc) = 1.286(13) GeV. (34) 

Our result agrees within the uncertainties with our last determination |10j|, but is consid- 
erably more precise. Using the three-loop relation [HI |42l |43] between pole- and MS-mass 
this corresponds to M^^~^°"^^ = 1.666 GeV. We refrain from providing an uncertainty for 
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1441 


1.224 ± 0.017 ± 0.054 


fit to B decay data, al(3o 




1.29 ± 0.07 


NNLO moments 


[36] 


1.319 ± 0.028 


lattice, quenched 


m 


1.26 ± 0.04±0.12 


lattice, quenched 


w\ 


1.301 ± 0.034 


lattice (ALPHA), quenched 




1.304 ± 0.027 


low-moment sum rules, NNLO 


[TO] 


1.25 ± 0.09 


PDG 2006 


[13J 



Table 8: Predictions for mc{mc). In the second column some keywords concerning the 
used method are given (NNLO: next-to-next-to-leading order; NNNLO: next-to-next-to- 
next-to-leading order) . 



n 


MO) 




Mil) 




/^(21) 




M30) 


r(3i) 


M32) 


M33) 


1 


1.0667 


2.5547 


2.1333 


3.1590 


3.4425 


0.0889 


-7.7624 


-0.0599 


1.5851 


-0.0543 


2 


0.4571 


1.1096 


1.8286 


3.2319 


5.0798 


1.9048 




4.0100 


7.2551 


0.1058 


3 


0.2709 


0.5194 


1.6254 


2.0677 


4.5815 


3.3185 




5.6496 


13.4967 


2.3967 


4 


0.1847 


0.2031 


1.4776 


1.2204 


3.4726 


4.4945 




3.9381 


17.2292 


6.2423 



Table 9: Coefficients of the photon polarization function in the MS scheme as defined in 
Eqs. (in]) and ([IS]) forn/ = 5. 



Mc since it is well known that the perturbative series between the MS and the pole mass 
is badly behaved. 

A comparison with a few selected determinations (published 2001 or later) is shown 
in Tab. [Hand Fig. [7] Within their respective errors all results are consistent. Let us also 
mention that a first investigation of the numerical effect of the four-loop results for the 
moments has been performed in Ref. [15] where basically the analysis of Ref. [10] has been 
adopted with updated parameters. The result of Ref. ^5] reads mc(mc) = 1.295(15) GeV. 

5 The bottom quark mass 

The same approach is also applicable to the determination of rrib. The coefficients (7„ are 
listed in Tab. [9l For the bounds on C^^^ we assume 

- 8.0 < °^ < 9.5 , -8.0 < °^ < 8.3 , -8.0 < °^ < 7.4 , (35) 

where the same criteria as in the charm case have been applied. The coefficients in Tab. [9] 
determine the theoretical moments through Eq. (l2Tl) where Qc has to be replaced by 
Qh = —1/3. The non-perturbative terms from the gluon condensate are small and will be 
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Figure 7: Comparison of recent determinations of mc{m^ (see also Tab. [8]). 





T(1S) 


T{2S) 


T(35) 


T{AS) 


MT(GeV) 
ree(keV) 
{a/a{Mr)f 


9.46030(26) 
1.340(18) 
0.932069 


10.02326(31) 
0.612(11) 
0.93099 


10.3552(5) 
0.443(8) 
0.930811 


10.5794(12) 
0.272(29) 
0.930093 



Table 10: Masses, electronic widths and effective electromagnetic couplings [IH] at s = M|- 
for the narrow T-resonances. 



neglected. We included the charm mass terms [19] of order (mc/m^)^ to the moments Cn 
which induce a shift of about —1 MeV in the bottom quark mass. 

The experimental results for the moments are listed in Tab. [TTJ The contribution from 
the resonances include T(IS') up to T(4S'). The values for the masses of the T resonances 
and their electronic width have been taken from Ref. |13| and are listed in Tab. [TOl The 
errors from the three lowest T resonances have been combined linearly, since the PDG 
values for the electronic widths are dominated by the measurement from CLEO [50J. 
The result is then combined in quadrature with the contribution from the T(45') reso- 
nance. The treatment of the bottom threshold region is similar to the one of the charm 
region. Measurements of R from threshold up to 11.24 GeV have been performed by 
the CLEO Collaboration more than 20 years ago [51j, with a systematic error of 6%. 



20 



n 




thresh 




A^exp 




Xl0(2-+1) 


Xl0(2"+1) 


Xl0(2n+1) 


Xl0(2n+1) 


1 


1.394(23) 


0.296(32) 


2.911(18) 


4.601(43) 


2 


1.459(23) 


0.249(27) 


1.173(11) 


2.881(37) 


3 


1.538(24) 


0.209(22) 


0.624(7) 


2.370(34) 


4 


1.630(25) 


0.175(19) 


0.372(5) 


2.178(32) 



Table 11: Moments for the bottom quark system. 



No radiative corrections were applied. The average value derived from the four data 
points below threshold amounts to ^ = 4.559 ± 0.034(stat.) which is 28% larger than 
the prediction from pQCD. However, a later result of CLEO [26j at practically the same 
energy, i?(10.52 GeV) = 3.56 ± 0.01 ± 0.07, is significantly more precise and in perfect 
agreement with theory. Applying a rescaling factor of 1/1.28 to the old CLEO data not 
only enforces agreement between old and new CLEO data and pQCD in the region below 
the T(4S'), it leads, in addition, also to excellent agreement between theory and experi- 
ment above threshold around 11.2 GeV where pQCD should be applicable also to bottom 
production. Further support to our approach is provided by the CLEO measurement of 
the cross section for bottom quark production at a/s = 10.865 GeV which is given by 
(rb{^/s = 10.865 GeV) = 0.301 ±0.002 ±0.039 nb [52]. The central value can be converted 
to i?fe(10.865 GeV) = 0.409. On the other hand, if one extracts i?b(10.865 GeV) from the 
rescaled CLEO data from 1984 ^ one obtains i?fe(10.865 GeV) = 0.425 which deviates 
by less than 4% from the recent result In Fig. [8] the original and the rescaled data 
from [51] is shown and compared to pQCD and data point from [26]. We thus extract the 
threshold contribution to the moments from the interval 10.62 GeV < y/s < 11.24 GeV 
by applying the rescaling factor to the data, subtract the "background" from u, d, s and 
c quarks and attribute a systematic error of 10% to the result. 

To evaluate the contributions to the moments from above 11.24 GeV, i.e. TW™"^', we 
use the prediction from pQCD for Rb{s). 

The results for the moments are listed in Tab. [TH those for mb(10 GeV) and mbirrib) 
in Tab. [121 Just as in the charm case, a remarkable consistency and stability is observed. 
For n = 1 the error is dominated by the experimental input. For n = 3 we obtain ±0.010 
from the experimental input, ±0.014 from and ±0.006 from the variation of fi. 

The sensitivity to the inclusion of higher orders is displayed in Fig. [91 For the lowest 
moment the error is dominated by the experimental uncertainty — nevertheless the theory 
error is reduced also in this case and the prediction stabilized. In general a significant 
improvement and stabilization is observed. 

The three results based on n = 1,2 and 3 are of comparable precision. The relative 
size of the contributions from the threshold and the continuum region decreases for the 
moments n = 2 and 3. On the other hand, the theory uncertainty estimated from the 
variation of /i and SC!^^^ is still acceptable. We therefore take the result from n = 2 (which 
is roughly between the n = 1 and n = 3 values) and add the uncertainty from "total" 
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Figure 8: In the upper plot the data for i?(s) are shown as published in Refs. |51| (circles) 
and |26] (triangle). The black curves are the predictions from pQCD outside the resonance 
region. In the lower plot the older data from [51] are rescaled by a factor 1/1.28. 

(cf. Tab. [T2|) and the one induced by ^C^^"-* linearly which leads to an error estimate of 
±25 MeV: 



7716(10 GeV) = 3.609(25) GeV, 



(36) 
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Figure 9: mf,(10 GeV) for n = 1, 2, 3 and 4. For each value of n the results from left to 
right correspond the inclusion of terms of order a], and in the coefficients Cn 
(cf. Eq. (JTSll ). Note, the central value for n = 4 and the uncertainties for n = 3 and 
n = 4 can not be determined in the case where only the two-loop corrections of order as 
are included into the coefficients Cn as the corresponding equation cannot be solved for 
mb(10 GeV). 



mb{mb) = 4.164(25) GeV. (37) 

This result can also be converted to a pole mass HH ia is] of Mf = 4.800 GeV. 
Again we do not display the additional error from the MS-pole-mass conversion, which is 
significantly larger. 

A comparison with a few selected determinations is shown in Tab. [13] and Fig. [TOl In 
Ref. [15] mbirrih) = 4.205(58) GeV was obtained on the bases of the four-loop result for 
n = 1 and the phenomenological analysis of Ref. [10] with updated resonance parameters. 

For various applications, either related to Z-boson decays or in connections to Grand 
Unified Theories (GUTs) the value of mb{fi) at high scales is of interest. The result for 
fi = Mz can be obtained in a straightforward way 

mb{Mz) = 2.834 ±0.019 ±0.017, (38) 

where the first error originates from the combined error listed in Eq. (!37l) . the second 
error from Sag- (We do not quote a theory error at this point.) 
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n 


mb(10 GeV) 


exp 


Us 


/i 


total 






1 


3.593 


0.020 


0.007 


0.002 


0.021 




4.149 


2 


3.609 


0.014 


0.012 


0.003 


0.019 


0.006 


4.164 


3 


3.618 


0.010 


0.014 


0.006 


0.019 


0.008 


4.173 


4 


3.631 


0.008 


0.015 


0.021 


0.027 


0.012 


4.185 



Table 12: Results for mb(10 GeV) in GeV obtained from Eq. (!23l) (adopted to the bottom 
quark case). The errors are from experiment, as and the variation of /x. The error from 
the yet unknown four-loop term is kept separate. Last column: central values are shown 
for mb{mb). 



mb{mb) (GeV) 


Method 


Reference 


4.164 ± 0.025 


low-moment sum rules, NNNLO 


this work 


4.19 ± 0.06 


T sum rules, NNLL (not complete) 


m 


4.347 ± 0.048 


lattice (ALPHA), quenched 




4.20 ± 0.04 


fit to B decay distribution, a'^Po 


m 


4.25 ± 0.02 ± 0.11 


lattice (UKQCD) 


m 


4.33 ± 0.10 


lattice, quenched 


m 


4.346 ± 0.070 


T(l^), NNNLO 




4.210 ± 0.090 ± 0.025 


T(15), NNLO 


[56] 


4.191 ± 0.051 


low-moment sum rules, NNLO 


[in] 


4.17 ± 0.05 


T sum rules, NNLO 


m 


4.20 ± 0.07 


PDG 


m 



Table 13: Predictions for mb{mb). In the second column some keywords concerning the 
used method are given (NNLO: next-to-next-to-leading order; NNNLO: next-to-next-to- 
next-to-leading order; NNLL: next-to-next-to- leading- logarithmic order). 



In connection with Yukawa unification in Grand Unified Theories the value mb{fi = rrit) 
in the nj = 6 theory may be of interest. Given Mt = 171.4 ± 2.1 [57], which we identify 
within the present uncertainties with the pole mass, we find 

mt{mt) = 161.8 ±2.0 ±0.2, (39) 

to three-loop accuracy [3T[ |121 US] . (Again we do not specify a theory uncertainty. The 
difference between the two- and three- loop result amounts to 426 Me V. ) Evolving rrib from 
/i = 10 GeV to mt{mt) and matching to the Uf = 6 theory we find 

mfe(161.8) = 2.703 ±0.018 ±0.019, (40) 



where the first error reflects the combined error from Eq. fl37l) and the second one the 
uncertainty due to a^. As stated above, the ratio 

\ [ = 59.8 ±1.3, 41 

mb{mt) 
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III 


1 1 1 1 1 1 1 1 1 1 1 1 1 
Kuehn, Steinhauser, Sturm 07 






Pineda, Signer 06 


' 1 




Delia Morte et al. 06 


^ — ' 


1 ■ 1 


Buchmueller, Flaecher 05 




1 


Mc Neile, Michael, Thompson 04 


^ ' — 




1 deDivitiis et al. 03 




■ 


1 Penin, Steinhauser 02 


' ^ ■ 


■ 

— 1 


Pineda 01 






Kuehn, Steinhauser 01 






Hoang 00 


1 ' ' 




PDG 2006 



4.1 4.2 4.3 4.4 4.5 4.6 4.7 



mb(m,) 

Figure 10: Comparison of recent determinations of rribirn^y) (see also Tab. [T3|l . 
should be a useful input for Grand Unified Theories. 

6 Conclusions 

A new determination of the charm and bottom quark masses has been presented. It is 
based on new experimental information on the electronic width of the narrow quarkonium 
resonances and on the cross section for charm and bottom production in the respective 
threshold region. In addition the analysis profits from the improved determination of 
from a new analysis of recent experiments. Four-loop results for the moments lead to a 
significant reduction of the theory error. The result can be directly expressed in terms 
of running MS masses conveniently expressed at a scale of 3 GeV and 10 GeV, respec- 
tively. Our final resuhs mc(3 GeV) = 0.986(13) GeV and mfe(10 GeV) = 3.609(25) GeV, 
correspond to mdrnc) = 1.286(13) GeV and mb{mf,) = 4.164(25) GeV. These values are 
consistent with our previous determination [10] but considerably more precise. Let us 
stress that, since rric is relatively small, mc{mc) shows a less stable behaviour against 
higher order corrections. Thus we propose to use mc(3 GeV) as the better choice for 
comparisons of the various methods. The evaluation of the running bottom mass at the 
scale of nit leads to mt{mt) /mf,{mt) = 59.8 ± 1.3, a result of interest for theories with 
Yukawa unification. We want to mention that for the charm quark mass our analysis 
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constitutes the only one where NNNLO corrections from theory are incorporated. In the 

case of the bottom quark there is only one further NNNLO analysis which is based on 
the mass of the T(IS') system. However, the accuracy of this approach is strongly limited 
due to large non-perturbative contributions. 
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